Anderson localization of polaron states 
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Using the vanishing of the typical polaron tunneling rate as an indicator of the breakdown of 
itinerancy, we study the localization of polaron states in a generic model for a disordered polaronic 
material. We find that extremely small disorder causes an Anderson localization of small polaron 
states. However, the ratio between the critical disorder strength needed to localize all states in the 
polaron band and the renormalized bandwidth is not necessarily smaller than for a bare electron. 

PACS numbers: 71.38.-k, 72.10.Di, 71.35.Aa 



The polaron concept has played a central role in the 
analysis of strongly coupled electron-phonon (EP) sys- 
tems since the pioneering work of Landau. Neverthe- 
less, even fundamental properties of a polaron are not 
yet completely understood, especially with respect to 
the effect of disorder. Polaron signatures are currently 
discussed in compounds such as the colossal magnetore- 
sistance manganites |2], the charge ordered nickelates 
and bismuthates and molecular organic semiconduc- 
tors -all complex materials, where imperfections po- 
tentially play an important role. 

The essence of polaron formation is the crossover from 
a weakly (phonon) dressed charge carrier to a strongly 
mass-enhanced, less mobile polaronic quasiparticle with 
increasing EP coupling strength. [|| |^ This process 
is usually called self-trapping, since the distortion that 
tends to bind the carrier in a deformable lattice is in- 
duced by the carrier itself, i.e., the trapping potential 
depends on the state of the carrier. [Q Nonetheless, even 
if the carrier is "confined" to a single lattice site (small 
polaron), this type of trapping does not imply localiza- 
tion. The tunneling between different lattice sites is still 
relevant and a self-trapped carrier resides in an itinerant 
polaron state. On the other hand, it is well known 
that for sufficiently strong randomness the itinerancy of 
(single electron) states is destroyed and the states are 
Anderson localized. ^ Obviously, disorder will affect 
the itinerancy of polaron states. Likewise, it was sug- 
gested that the "mobility edge" , separating itinerant (ex- 
tended) and localized states, will be strongly influenced 
by the EP interaction At the moment, however, 

despite its apparent relevance for various transport pro- 
cesses in polaronic materials, little is known rigorously 
about how Anderson localization affects polaron states, 
and vice versa. 

As a first step toward addressing this problem, in 
Ref. ijlll a single polaron in a Holstein model with site- 
diagonal, binary-alloy-type disorder was studied within 
the dynamical mean field approximation (DMFA). 
The DMFA, however, cannot (fully) discriminate be- 
tween itinerant and localized states (as shown below), 
mainly because the randomness is treated at the level of 
the coherent potential approximation. |T3] 



The purpose of this paper is to remedy this shortcom- 
ing, and to present an advanced theory that is capable 
of describing both the Anderson localization and self- 
trapping phenomena. As a generic model we consider 
the single-particle Anderson-Holstein Hamiltonian 

n = ^e^n, - J^(4cj +H.c.) 

-^/E;nY,ib^ + bl)n, + nY,blh, (1) 

i i 

where J is the electron transfer amplitude between neigh- 
boring sites, fl denotes the bare frequency of the optical 
phonon (h — 1), and Ep is the polaron shift. In Eq. (|^), 
the on-site energies {e^} are assumed to be independent, 
identically distributed random variables with a probabil- 
ity density p{e,) = (1/7)^(7/2 - |e,|). 

Without disorder, the physical properties of the Hol- 
stein model are governed by three parameter ratios: the 
adiabaticity a = fl/ J [determining whether the electron 
is light (a ^ 1) or heavy (a » 1) (Ref. 0)], and two 
dimensionless EP coupling constants A = Ep/2J and 
= Ep/il. Extensive numerical studies of the Hol- 
stein model on finite one-dimensional (ID) and 2D lat- 
tices have shown that polaron formation sets in if two 
conditions A > 1 and g^ > 1 are fulfilled. H The 
internal structure of the polaron depends on A and a. As 
a result, disorder affects polaron states quite differently 
in the adiabatic (a <C 1), nonadiabatic (a 1), and an- 
tiadiabatic (a ^ 1) cases. Furthermore, in the crossover 
region, where polaron formation sets in, long-range tun- 
neling induced by the EP coupling ||l5| is expected to 
decrease the tendency toward localization. 

To analyze the effects of disorder on polaron forma- 
tion quantitatively, we adopt the recently developed sta- 
tistical DMFA (statDMFA) (Ref. 0) to the Anderson- 
Holstein model. This method is an extension of the self- 
consistent theory of localization of bare electrons to 
models with interactions, which has been successfully ap- 
plied, e .g.^ to the disordered Hubbard and Anderson 
lattice 1 19 models. The statDMFA is essentially a proba- 
bilistic method, based on the self-consistent construction 
of random samples for the (local) physical quantities of 
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interest. 

The statDMFA is conveniently constructed for a Bethe 
lattice with finite connectivity K . For later purposes, we 
rescale the transfer amplitude J J j^fK and define A = 
E^I2J and a = ^13. Applying the statDMFA (Ref. |l7|) 
to model (|l|), the local zero-temperature single-electron 
Green function reads Qi{z) = [z — — Hi(z) — T,i{z)]^^ 
{z = (jj + irj), with the EP self-energy and the 

hybridization function 



age. In that case, the geometric mean (typical value) 



K 



i/,(z)-E,(z)- 
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In the spirit of the self-consistent theory of localiza- 
tion, the statDMFA ignores that the functions on 
the right-hand side of Eq. (||) should be calculated for 
the Bethe lattice with the site i removed. It furthermore 
takes K as the typical number of terms, although for the 
initial site there would be K+l terms. As a consequence, 
the bare local density of states (LDOS) becomes semiel- 
liptic (with a half-band- width Wq/2 = 2J), irrespective 
of the value of K, which makes the mean field character 
of the StatDMFA evident. |20) 

The most severe approximation within the statDMFA 
is related to the calculation of Si(z). Although the con- 
nectivity of the Bethe lattice is finite in Eq. (g) , the self- 
energy is determined as if the Bethe lattice had infinite 
connectivity. |^ The self-energy is then local (as in the 
DMFA), |n|, H and reads S,(z) = Epin/{[Fl^\z)]-^ - 
Ep2n/...) with [F^P\z)]-^ = z-pn~ h[^\z) and 
h[^\z) ^ Hi{z — pft). Here the energy shift keeps track 
of the number of virtual phonons, and the index p is a 
positive integer, because of the absence of thermally ex- 
cited phonons. 

A natural measure of the itinerancy of a polaron state 
is the tunneling rate from a given site defined as the 
imaginary part of the hybridization function Ti{uj) — 
{TTjyK)J2f^,Nj{uj), where N,{uj) = -(l/^)Im5, (^) 
is the LDOS. A vanishing tunneling rate Ti{u}) implies a 
localized state at energy lo. 

Due to the randomness in the on-site energies, the tun- 
neling rate Ti{Lu) is a random variable, and the ques- 
tion of whether it vanishes or not depends on the prob- 
ability density exhibiting, quite generally, characteris- 
tic differences for itinerant and localized states. 
For an itinerant state, the probability density is narrow 
and (approximately) symmetric. Thus, the most proba- 
ble value and the arithmetic mean (average), r'*^(u;) = 
{'i^/N)^f Ti{u}), N being the total number of lattice 
sites, are almost identical and r^^(w) will be a good es- 
timate for the tunneling rate. On the other hand, for a 
localized state, the probability density is very wide with 
an extremely long tail toward large values, implying that 
the most probable value is much smaller than the aver- 



r*yP(tj) = exp 
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seems to be a more appropriate estimate of the tunnel- 
ing rate. Accordingly, a vanishing typical (rather than 
average) tunneling rate should be taken as a localization 
criterion, jl^ 

The typical tunneling rate is most directly calculated 
from a recursion scheme for the hybridization function 
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Epin 



(4) 



Ep2^ 



Demanding that Hj^^\z) are independent, identically dis- 
tributed random variables, Eq. (|^) can be used to con- 
struct a random sample {H^^\z)} (i = l,...,iV and 
p — 0,1, ...,M) for the hybridization function, starting 
from an initial random configuration which is successively 
updated with a sampling technique similar to the one de- 
scribed in Ref. ||l^. More precisely, we are constructing 
a self-consistent random sample because the probability 
density for H'f\z) is assumed to be the same on both 
sides of Eq. (^. Clearly, the larger N the better the 
statistics. A correct description of the polaron effect re- 
quires the cutoff M for the continued fraction to be much 
larger than the average number of phonons bound in the 
polaron. Specifically, for the results presented in this pa- 
per = 50 000 and M = 35. 

It is instructive to consider Eq. (||) for K ~ oo. Apply- 
ing the law of large numbers shows that H'f^ (z) can be 

replaced by its average H''''^p\z). Then, (1/ K)J2f (...) 
denotes the (sample) average solely over e;, and Eq. (^) 
reduces to the DMFA equation for H'^^^p^z). (The scal- 
ing of J makes this point particularly clear.) Obviously, 
the tunneling rate is then 7tJ'^N'^^{z), i.e., it is propor- 
tional to the average LDOS, which is always finite, ir- 
respective of whether the state at energy lo is itinerant 
or locahzed. §, § Thus, for K = oo, i.e., in the DMFA, 
itinerant and localized states cannot be distinguished. 

In the following we present representative numerical 
results for K = 2 and rj = 10~^. We start our discussion 
with the (nonadiabatic-to-antiadiabatic) weak EP cou- 
pling regime. Figure |l| displays typical tunneling rates 
r^^P{uj;j) as a function of disorder strength 7 for two 
energies lo, one below and one above the optical-phonon 
emission threshold, which, due to the large value of a, 
lies far inside the band (see the inset) . Below the thresh- 
old, Iml](cj) — and the states are coherent. Above the 
threshold, however, inelastic polaron-phonon scattering 
yields ImI](Lj) =/= 0, i.e., the states are incoherent, im- 
plying that the quantum interference needed for localiza- 
tion is significantly suppressed. We therefore expect 
the critical disorder strength 7c, for which T^'^^{lo; 7) van- 
ishes, to be larger above the phonon emission threshold 
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FIG. 1: Typical tunneling rates r*^P(o;;7) in the (non-to- 
antiadiabatic) weak EP coupling regime (Wo ~ 1-0, a = 1.2, 
and A = 0.067). Filled circles denote data for A = 0. The 
inset shows the LDOS N{ui) and the imaginary part of the 
EP self-energy ImE(aj) for 7 = 0; dashed and (vertical) dotted 
lines indicate, respectively, the density of states for A = and 
the energies ui for which r'^P(tj;7) is plotted. 
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FIG. 3: The upper and lower panels show, respectively, typi- 
cal tunneling rates T^^'^ (lu;')) in the (adiabatic) intermediate 
(3 = 0.4, A = 1.8) and (non-to-antiadiabatic) strong EP cou- 
pling regime (q = 1.2 and A — 6.66) using Wo = 2.0. The 
insets display the subband LDOS A*'(a;) and the imaginary 
part of the EP self-energy ImE (a;) for 7 = (vertical dotted 
lines indicate the energies lo for which T^^^{lu;^) is plotted). 
Note the small spectral weight of the lowest polaron subband. 



FIG. 2: Evolution of the LDOS N{lj) (thick line) and the 
imaginary part of the EP self-energy ImE(tj) (thin line) for 
the ordered Holstein model (7 = 0) with increasing EP cou- 
pling A = 0.4,0.7, and 1.0 (from top to bottom). Results are 
given for a — 0.4 and Wo = 1.0. 



than below it. Indeed, we find (7c/Wo)w=-o.4 ~ 2.25 
and (7c/W^o)w=o.4 ~ 3.0. It is rather encouraging that 
the statDMFA recovers this basic feature of the physics 
of localization. Since the EP coupling is weak the states 
below the threshold are basically bare electron states and 

(7c/W^o)^=-o.4 ~ {lc/Wo)t=\A ~ 2.0, the correspond- 
ing ratio for a bare electron. 

Increasing the EP coupling, polaron formation starts 
and the LDOS fragments into an increasing number of 
(more or less) coherent subbands. Figure || illus- 
trates, for the ordered Holstein model (7 = 0), the evo- 



lution of the LDOS with increasing EP coupling A for 
fixed a and Wq- As far as localization is concerned, the 
lowest subband is of particular interest, because it is com- 
pletely coherent and the quantum interference needed for 
localization is not masked by inelastic polaron-phonon 
scattering. Numerical studies of the ordered Holstein 
model demonstrated that the band dispersion of the low- 
est subband depends on both A and a, and may differ 
significantly from the rescaled bare band, jlj, |l^ Specif- 
ically, the states at the bottom of the subband are mostly 
electronic and rather mobile due to long-range tunneling 
induced by EP coupling, whereas the states at the top 
of the subband are rather phononic and immobile. |15| ] 
Therefore, we expect that the critical disorder strength 
needed to localize states at the top of the subband is 
much smaller than for states at the bottom. This behav- 
ior should be most pronounced in the (adiabatic) inter- 
mediate EP coupling regime. 
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The upper panel of Fig. showing typical tunneling 
rates r*yp(a;;7) as a function of disorder strength 7 for 
the lowest polaron subband, impressively confirms this 
striking behavior. The strong asymmetry of the LDOS 
(see the inset) is a direct consequence of the "hybridiza- 
tion" with the (optical) phonon branch, leading to a band 
flattening at the Brillouin-zone boundary. As can be 
seen, states at the zone boundary are very susceptible to 
disorder and the critical disorder strength 7c for which 
r*^P(a;;7) vanishes is substantially smaller than at the 
bottom of the subband. Clearly, the absolute scale of 
disorder affecting the subband is on the order of the sub- 
band width Wi, which is strongly renormalized in com- 
parison with the bare band width Wo- Different from 
the pure Anderson model, the critical disorder strength 
needed to localize all states of the subband is determined 
by the states at the bottom (not at the center). We 
find (7c /W^i) complete ~ 2.8, which is larger than the cor- 
responding ratio for a bare electron [(7c/VFo)compictc ~ 
2.25] because of the long-range tunneling induced by the 
EP coupling. Thus, contrary to naive expectations, in 
the intermediate-coupling (crossover) regime, a polaron 
is more difficult to localize than a bare electron. 

Finally, in the lower panel of Fig. we present typ- 
ical tunneling rates r*yp(u;;7) in the (nonadiabatic-to- 
antiadiabatic) strong EP coupling regime. The phonon 
admixture of the polaron states within the lowest sub- 
band is now nearly energy independent. Concomitantly, 
long-range tunneling as well as band flattening are neg- 
ligible and the LDOS is rather symmetric (see the in- 



set). [|T^ Accordingly, disorder should affect all the 
states within the subband more or less equally. In- 
deed, the critical disorder strength 7c for which r^^P{uj; 7) 
vanishes is now almost energy independent. We find 
(Tc/W^i) complete ~ 2.5, which is smaller than the ratio for 
intermediate EP coupling, indicating that in the extreme 
strong EP coupling adiabatic (A ^ 1) and nonadiabatic- 
to-antiadiabatic {g^ ^ 1) regimes (7c/W^i)compiete ap- 
proaches the corresponding ratio for a bare electron. 
Here, the band collapse does only change the overall en- 
ergy scales. 

To summarize, we presented, as a first step toward 
a transport theory for disordered polaronic materials, a 
theoretical approach capable of describing polaron for- 
mation and Anderson localization. We showed that the 
typical tunneling rate from a given site is an appropriate 
quantity in order to distinguish itinerant and localized 
polaron states. Focusing on the disorder-induced local- 
ization of a single polaron in the intermediate-to-strong 
EP coupling regime, we found that extremely small ma- 
terial imperfections turn itinerant polaron states into lo- 
calized polaron states, suggesting that in real materials 
(small) polarons are presumably always localized. In the 
crossover regime, the critical disorder strength, in units 
of the subband width, needed to localize all states of the 
lowest polaron subband can be larger than that for the 
pure Anderson model. From this perspective, disorder- 
induced localization of a polaron is not necessarily easier 
than localization of a bare electron. 
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